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ABSTRACT 


Amplitude  evolution  equations  are  derived  for  viscous  gravity  waves  and  for 
viscous  capillary-gravity  waves  with  surfactants  in  water  of  infinite  depth.  Multiple 
scales  are  used  to  describe  the  slow  modulation  of  a  wave  packet,  and  matched 
asymptotic  expansions  are  introduced  to  represent  the  viscous  boundary  layer  at 
the  free  surface.  The  resulting  dissipative  nonlinear  Schrodinger  equations  show 
that  the  largest  terms  in  the  damping  coefficients  are  unaltered  from  previous 
linear  results  up  to  third  order  in  the  amplitude  expansions.  The  modulationad 
instability  of  infinite  wavetrains  of  small  but  finite  amplitude  is  studied  analytically 
and  computationally.  For  capillary-gravity  waves  a  band  of  Weber  numbers  is 
found  in  which  the  linear  analysis  guarantees  neutral  stability  in  the  absence  of 
viscous  dissipation.  The  corresponding  spectral  computation  shows  modulation 
features  that  represent  a  sm2Jl- amplitude  recurrence  not  directly  related  to  the 
Benjamin-Feir  instability. 
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1  Introduction 


Although  water  waves  are  adequately  described  by  inviscid-flow  theory  in  most  cases,  the 
effects  of  viscosity  and  interfacial  properties  become  important  for  waves  of  short  wavelength. 
For  sufficiently  short  waves,  the  shear-stress  boundary  condition,  which  is  neglected  in  classical 
studies,  must  then  be  satisfied  at  the  free  surface.  If  the  free  surface  is  clean  and  the  air  above 
is  ignored,  the  shear  stress  vanishes.  However,  when  a  layer  of  contaminant  (surfactant)  is 
present  on  the  free  surface,  its  concentration  varies  in  accordance  with  the  motion  of  the  free 
surface,  causing  a  surface-tension  gradient  that  must  be  balanced  by  a  nonzero  surface  shear 
stress.  Therefore,  the  characteristics  of  the  surfactant  must  be  examined  before  the  boundary 
conditions  can  be  prescribed  completely. 

The  highly  dissipative  effect  of  surfactants  has  been  observed  since  classical  times  (see  Pliny, 
77  A.D.,  for  example).  A  qualitative  explanation  of  its  cause  in  terms  of  the  variations  in 
surface  tension  was  given  by  Reynolds  (1880),  as  cited  by  Lamb  (1932).  Levich  (1941,  1962) 
later  performed  a  more  thorough  analysis,  and  an  extension  of  this  work  was  described  in  the 
review  of  Lucassen-Reynders  ic  Lucassen  (1969).  In  these  mathematical  models  the  surfactant 
was  considered  as  an  interfacial  film  of  infinitesimal  thickness.  Closed-form  solutions  for  the 
damping  coefficient  were  given  for  linear  capillary  waves,  suid  the  damping  effect  was  shown 
to  increase  dramatically  from  that  for  a  clesui  surface.  A  series  of  more  complicated  models 
proposed  by  Goodrich  (1960,  1962)  treated  the  surfactant  as  an  anisotropic  rheological  body.  A 
simpler  viscoelastic  interfacial  constitutive  equation  that  compared  favorably  with  experimental 
measurements  was  developed  by  Addison  &  Schechter  (1979).  Extensive  references  describing  the 
various  effects  of  surfactants  on  water  waves  were  given  by  Lucassen  (1982)  and  Garrett  (1986). 

In  the  aforementioned  investigations,  the  study  of  viscous  damping  associated  with  water 
waves  has  been  confined  to  linear  theory,  in  which  the  wave  amplitude  is  infinitesimal  and  the 
free-surface  boundary  conditions  are  applied  at  a  mean  surface.  In  the  present  work,  we  focus 
on  the  wave  nonlinearity  and  adopt  the  simple  elastic  model  of  Lucassen-Reynders  Sc  Lucassen 
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(1969)  and  Lucassen  (1982),  which  is  suitable  for  many  insoluble  surfactants. 

The  usual  theory  of  weakly  nonlineair  waves  is  based  on  the  assumptions  that  the  fluid  is 
inviscid  and  the  motion  is  irrotational.  Of  particular  interest  are  solutions  that  address  the 
evolution  of  a  wave  packet  and  the  instability  of  a  Stokes  wavetrain.  Benjamin  &  Feir  (1967) 
and  Whitham  (1967)  demonstrated  in  different  ways  that  progressive  waves  of  finite  amplitude 
on  deep  water  are  unstable.  The  evolution  equation  for  the  slowly  varying  amplitude  is  a  non¬ 
linear  Schrodinger  equation  (NLS),  first  introduced  for  water  waves  by  Zakharov  (1968)  with 
a  variational  method  and  obtained  by  Hasimoto  ti  Ono  (1972)  using  the  method  of  multiple 
scales.  Peregrine  (1983)  summarized  various  analytical  solutions  to  the  NLS  equation,  and  Yuen 
ii  Ferguson  (1978),  Weideman  k.  Herbst  (1987),  and  many  others  studied  numericaJ  solutions. 
An  extension  to  three-dimensional  wave  packets  was  given  by  Davey  k  Stewartson  (1974),  and 
capillary  effects  were  considered  by  Djordjevic  k  Redekopp  (1977).  The  damping  effect  has  been 
included  in  the  two-dimensional  NLS  equation,  e.g.,  by  Pereira  k  Stenflo  (1977),  through  the 
addition  of  a  conjectured  small  dissipation  term.  A  primary  objective  of  the  present  work  is 
to  develop  a  systematic  derivation  of  the  amplitude  equation  in  the  presence  of  viscosity  and 
a  thin  surfactant  film.  The  result  is  a  dissip2d;ive  nonlinear  Schrodinger  equation  (dissipative 
NLS),  which  contains  a  linear  dissipation  term;  the  other  terms  remain  unaltered  from  their 
counterparts  for  inviscid  flow. 

We  adopt  the  multiple-scale  method  to  represent  the  wave  modulation  and  a  boundary-layer 
scaling  to  incorporate  the  effects  of  viscosity  and  surfactant.  For  a  water  depth  that  is  large 
in  comparison  with  the  length  of  the  wave  packet,  a  slow  modulation  in  the  vertical  direction 
must  be  added,  as  introduced  by  McGoldrick  (1970)  and  explained  further  by  Mei  (1983).  The 
free-surface  boundar>  layer  has  been  considered  in  several  studies  of  the  mean  drift  in  water 
waves,  each  in  terms  of  coordinates  such  that  a  boundary-layer  variable  could  be  measured  from 
the  actual  location  of  the  free  surface.  Longuet-Higgins  (1953)  and  Liu  k  Davis  (1977)  used  a 
local  curvilinear  coordinate  system  moving  with  the  free  surface,  Phillips  (1966)  and  Craik  (1982) 
chose  a  conformal  transformation,  and  fjnluata  k  Mei  (1970)  introduced  a  Lagrangian  system. 
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In  the  present  study,  our  goal  is  somewhat  different,  and  the  analysis  is  extended  to  higher  order. 
One  choice  of  coordinates  that  we  considered  initially  was  an  orthogonal  system  obtained  from  a 
higher-order  conformal  transformation.  It  appeared  to  us,  however,  that  a  more  strsughtforward 
formulation  can  be  achieved  by  a  nonorthogonal  system  with  one  coordinate  measured  horizon¬ 
tally  and  the  other  measured  vertically  from  the  exact  location  of  the  free  surface. 

The  formulation  for  two-dimensional  capillary-gravity  waves  is  developed  in  §2,  where  the 
coordinate  system  and  governing  equations  are  presented.  Suitable  variables  for  the  methods 
of  multiple  scales  and  matched  asymptotic  expansions  are  introduced,  and  the  first-order  inner 
and  outer  solutions  are  derived.  In  §3,  a  dissipative  NLS  equation  is  derived  for  viscous  gravity 
waves  with  slowly  varying  amplitude.  A  solution  for  infinite  wavetrains  (the  viscous  equivalent 
of  Stokes  waves)  is  obtained,  and  the  resulting  damping  coefficient  is  determined.  In  §4,  capillary 
effects  are  added.  The  resulting  dissipative  NLS  equation  describes  the  effect  of  surfactants  on 
the  evolution  of  weakly  nonlinear  waves.  Nonuniformities  associated  with  the  second-  and  the 
third-harmonic  resonances  and  with  pure  capillary  waves  are  discussed.  Effects  of  viscosity  and 
surfactants  on  the  stability  of  weakly  nonline2ir  waves  are  analyzed  in  §5.  Numerical  solutions  to 
the  dissipative  NLS  equations  are  obtained  in  §6;  these  solutions  show  the  large-time  modulation 
of  viscous  gravity  and  capillary-gravity  waves  in  the  presence  of  surfactants. 

2  Formulation 

We  consider  the  slow  evolution  of  two-dimensional  surface  waves  influenced  by  an  insoluble 
surfactant  as  well  as  by  dispersion  and  weak  nonlinearity.  The  extent  of  the  wave  packet  is  taken 
to  be  small  in  comparison  with  the  fluid  depth,  so  that  the  depth  may  be  regarded  as  infinite. 
The  fluid  below  the  free  surface  has  constant  density  p  and  viscosity  p,  and  the  air  above  the  free 
surface  will  be  ignored.  The  surfactant  is  characterized  by  a  surface  dilational  modulus  M,  which 
measures  the  resistance  to  the  compression/expansion  type  of  surface  deformation  (Lucassen, 
1982).  Changes  in  surfactant  concentration  cause  the  surface  tension  <r  to  vary  in  time  and 
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space.  The  flow  is  caused  by  initial  surface  disturbances  with  small  but  finite  amplitude,  and  is 
irrotational  everywhere  except  in  the  boundary  layer  beneath  the  free  surfjw:e,  since  the  viscosity 
is  considered  to  be  small. 

In  general,  the  boundary-layer  thickness  may  be  much  smaller  than  the  wave  amplitude. 
Since  the  boundary-layer  part  of  the  solution  is  expressed  in  terms  of  an  inner  variable  measured 
from  the  surface,  matching  with  an  outer  potential-flow  solution  is  facilitated  by  the  introduction 
of  a  coordinate  system  that  moves  with  the  free  surface.  One  such  system  is  a  nonorthogonal 
coordinate  system  {x,y)  generated  simply  by  subtracting  the  free-surface  elevation  rj  from  the 
vertical  coordinate  of  a  Cartesian  system  (ar',  y'): 

x  =  x',  y  =  y' -*?(*.  0.  (1) 

where  Ti{x,t)  =  77(1',  <)  and  y  is  directed  upward  from  the  mean  free  surface.  The  two  components 
of  the  momentum  equation  are  obtained  by  applying  (1)  to  the  nondimensional  Navier-Stokes 
equations  expressed  in  the  Cartesian  coordinate  system  (x',y'): 

ut  +  uut  -h(v~rj,-  UTfg)uy  =  -pg  +  rjgpy  •+  e  [u,,  +  (!-!•  -  2r)gUgy  -  r;„u,]  (2) 

V,  +  UVg  +  {v  -  t},  -  Uf}g)Vy  =  -Py  +  f  (v„  -b  (1  -I-  Tlg)Vyy  -  2rjgVgy  “  TJggVy]  .  (3) 

Here,  the  reciprocals  of  the  wavenumber  k  and  the  frequency  u,  of  the  linear  fundamental  wave  are 
used  as  the  length  and  time  scales,  respectively.  The  horizontal  and  vertical  velocity  components 
u  and  V  are  scaled  by  Uo/k,  and  the  pressure  p  (with  the  hydrostatic  component  subtracted)  is 
scaled  by  pul/k"*.  Since  the  velocity  vector  is  still  written  using  the  base  vectors  of  the  Cartesian 
coordinate  system  (*',y'),  only  the  chain  rule  has  been  used  to  derive  (2)-(3).  The  reciprocal  of 
the  Reynolds  number,  which  is  the  primary  small  parameter,  is  defined  as 


iFYom  the  continuity  equation  for  constant  density,  we  have 

-b  Vj  =  0.  (4) 
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If  the  elements  of  the  nondimensional  stress  tensor  referred  to  the  {x' ,y')  coordinate  system 
are  denoted  by  Tij,  and  nj  and  nj  are  the  x'  and  y'  components  of  the  unit  vector  normal  to  the 
surface,  the  surface  force  per  unit  area  has  x'  and  y'  components  Fi  and  F2  given  by  Fi  =  TijUj. 
Since  the  slope  of  the  free  surface  is  t)*,  the  components  of  the  outward  normal  are  ni  =  — tJi/L 
and  no  =  l/L,  where 

L  =  {l  +  r,iy/\ 


The  normal  and  tangential  stresses  then  become  and  FiU,  where  ti  =  nj  zuid  =  —ni- 
The  normal  stress  on  the  free  surface  should  balance  the  product  of  the  curvature  eind  the 
surface-tension  coefficient  when  the  pressure  above  the  free  surface  is  taken  to  be  zero.  In  the 
present  coordinate  system,  this  condition  is  expressed  as 

2(  T 

P  = -[2  [vli^x  -  VxrJiy)  -  +  (I  +  T]l)vy] Cri  -  at  y  =  0,  (5) 
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where 


J-  /W.J* 


are  the  reciprocal  of  the  Rroude  number  and  the  Weber  number,  respectively.  Here,  g  is  the 
gravitational  acceleration  acting  vertically  downward,  and  is  the  undisturbed  value  of  the 
surface  tension.  The  boundary  condition  (5)  is  applied  at  the  exact  location  of  the  free  surface 
y  =  0  due  to  the  transformation  (1). 

The  tangential  stress  on  the  free  surface  balances  the  surface-tension  gradient  <t,,  which  can 
be  related  to  the  surface  strain.  If  the  displacement  of  a  material  point  on  the  surface  has  x'  aind 
y'  components  C  and  r/,  the  displacement  component  in  a  direction  along  the  surface  is  ^ti  +  17I3. 
The  surface  strain  equals  the  derivative  of  the  displacement  with  respect  to  distance  s  along  the 
surface,  and  the  gradient  of  surface  tension  is  assumed  to  be  proportional  to  the  derivative  of  the 
strain.  This  gives  a  nonlinear  extension  of  the  insoluble-surfactant  model  used  by  Levich  (1962), 
Lucassen-Reynders  k  Lucassen  (1969),  and  others; 

—  (C  +  nvx),. , 


where  the  surface  dilational  modulus  M  (or  inverse  surface  compressibility)  measures  the  change 
in  surface  tension  resulting  from  a  unit  fractional  area  change.  A  complete  rheological  interpre¬ 
tation  of  this  quantity  was  given  by  Lucassen-Reynders  ic  Lucassen  (1969).  The  nondimensional 
parameter  k^M/ipu^)  is  small  and  is  chosen  to  be  0(c^^^)  so  that  the  effect  of  surfactants 
will  appear  in  the  first  approximation  to  the  boundary-layer  solution  (later  in  this  section).  A 
nondimensional  surface  dilational  modulus  is  defined  as 

k^M 

In  the  present  coordinate  system,  the  required  surface  shear  stress  is  then 

2»j,  (uy  -  u,  +  fj,fUy) -I- (1  -  rjJ)(uy -I- V, -»?*t),)  =  €“‘/®/c(C-(- »?*?*)„  at  y  =  0.  (6) 

Since  both  the  horizontal  and  the  vertical  displacements  of  the  free  surface  are  introduced, 
two  kinematic  conditions  on  the  free  surface  are  needed  to  complete  the  governing  equations: 


t)  =  n*  +  utjg 

at  y  =  0 

(7) 

«  =  C«  +  uCx 

at  y  =  0. 

(8) 

Since  we  consider  weak  nonlinearity,  we  abo  require  the  initial  amplitude  a  of  the  fundamental 
mode  to  be  small.  The  relative  importance  of  nonlinearity,  surfactant,  and  viscosity  is  determined 
through  the  scaling  of  o,  k,  and  c,  respectively.  Here,  we  choose  k  to  be  0(1)  to  observe  the  effect 
of  surfactants  in  the  first  approximation.  The  relationship  between  a  and  the  boundetry-layer 
thickness,  which  is  0(e‘^^),  is  chosen  such  that  an  amplitude  equation  with  viscous  effects  can 
be  determined  at  third  order. 

We  introduce  the  slow  variables 

T  =  ah,  ^  =  o(*  -  Cft),  y  =  ay,  (9) 

where  Cy  denotes  the  group  velocity  of  the  primary  progressive  wave.  The  first  two  small  scales  are 
identical  to  those  used  by  Davey  &  Stewartson  (1974)  and  Djor^jevic  &  Redekopp  (1977)  for  an 
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inviscid  fluid  with  finite  depth.  The  slow  modulation  in  the  vertical  direction  is  added  to  suppress 
the  inconsistency  in  the  third-order  solutions  for  infinite  depth,  as  explained  by  Mei  (1983)  for 
an  inviscid  fluid.  The  slow  variation  in  the  horizontal  direction  thus  requires  a  corresponding 
long-scale  variation  in  the  vertical  direction,  and  the  nondimensional  depth  is  taken  to  be  much 
larger  than  the  extent  of  the  wave  packet,  which  is  0{a~*). 

The  flow  outside  the  boundary  layer  can  be  described  by  a  velocity  potential  which  should 
satisfy  the  Laplace  equation.  In  the  present  coordinate  system, 

~  (10) 

We  expand  the  velocity  potential  4  for  small  amplitude  a; 

4>  =  cl4i  +  47  +  (11) 

By  introducing  the  slow  variables  (9)  and  substituting  the  expansion  (11)  into  (10),  we  then 
obtain,  for  increasing  powers  of  a, 

=  0,  (12) 

<t>7tx  +  4>7n  —  (13) 

-*■  <^3yy  =  2*7ir(^2«y  +  +  ^l<y )  +  2(»72*  +  »7l()^lry  +  ^ly(»72r*  +  2»7ir{) 

+’7i«(^iy  +  ^Jy)  -  2^3x(  —  4iu  ~  2^2yy  -  4\fii  -  llx^iyy  (1^) 
The  lowest-order  equation  ( 12)  gives 

4i  =  {AE  +  A^E-^)e»,  (15) 

where 

£  =  exp{i(i  -  t)}  (16) 

and  A{^,  r)  describes  a  slow  variation  of  amplitude  with  a  mean  value  of  ^  at  r  =  0,  since  a  is 
the  initial  amplitude  of  this  mode.  The  superscript  c  denotes  the  complex  conjugate,  and  t  is 
Since  the  factor  e*  in  (15)  decays  exponentially  with  depth,  the  slowly  varying  amplitude 
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A  need  not  depend  on  y.  In  general,  a  long-scale  velocity  potential  ^i(^,y,r)  can  be  added  to 
the  right-hand-side  of  (15),  which  should  satisfy  -t-  =  0  to  avoid  secular  terms  arising 

from  (14).  Also,  the  boundary  condition  ^ij(y  =  0)  =  0  should  be  imposed  to  suppress  secular 
terms  in  the  second-order  free-surface  elevation  ijj  below.  In  the  present  study,  we  set  <t)°  =  0. 
An  example  of  nonzero  4>°  for  a  time-varying  current  is  given  by  Mei  (1983). 

The  velocity  components  are  obtained  from  the  gradient  of  the  velocity  potential  (11)  in  the 
x'  and  y'  coordinates  with  the  transformation  (1); 

u  =  ai{AE—A‘E  *)e*  4- a*(^2r  +  ~  r/ir^iy)  +  •  •  •  (17) 

V  =  a{AE  A^E  ^)e* -f  a* )-♦-■•■ .  (18) 

^From  the  Bernoulli  equation,  the  pressure  for  the  outer  region  of  ir rotational  flow  is 

p  =  ai(AE  —  A‘E~^)e^  —  a^[4>7t  —  —  Vu^iy  +  2^^**  ^i»)]  “  ~ 

+  (0?v  -  »72«  +  Cgrm)4>\y  +  (^5,  +  ~  +  •  •  •  •  (19) 


In  the  boundary  layer,  the  system  (2)-(8)  is  solved  after  introducing  the  slow  variables  (9) 
and  the  inner  coordinate 


y* 


(20) 


Velocity  components,  pressure,  and  free-surface  displacements  are  expanded  for  large  Reynolds 
number.  The  gauge  functions  for  the  inner  expansions  are  to  be  determined  after  the  scade  of  the 
amplitude  is  chosen.  Independent  of  the  scale,  the  first-order  boundary-layer  equations  become, 
in  the  order  of  usage; 


W,.  =0 

Tiu  =  vj  at  y*  =  0 

Pi,.  =  0 

pI  =  -  Trhgs  at  y*  =  0 

“i,*,*  ~  “I«  ~  pI* 


(21) 

(22) 

(23) 

(24) 

(25) 
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at  y*  =  0 


(26) 


«!»•«  =  «“l*r 

Cii  =  uj  at 


y*  =0, 


(27) 


where  uj,  and  pj  are  the  first-order  inner  solutions.  The  boundary  condition  for  uj  in  (26) 
is  obtained  by  differentiating  (6)  with  respect  to  time  and  applying  (8).  The  solutions  must  also 
match  with  the  outer  solutions  (17)-(19).  After  matching  conditions  are  applied,  the  solutions 
become 


=  AE  -1-  c.c. 

(28) 

rji  =  iAE  +  c.c 

(29) 

p*  =  iAE  +  c.c. 

(30) 

t»j  =  (Gj  i)AE  -f  c.c. 

(31) 

Cl  =  (»Go  —  l)AE  +  c.c.  , 

(32) 

where 


Qi  = 


v/2« 


1  -  »(1  -  V2«) 

and  c.c.  denotes  the  complex  conjugate.  The  linear  dispersion  equation 


G-l-r=  1 


(33) 


is  obtained  from  the  boundary  and  matching  conditions  for  the  pressure  and  is  not  affected  by 
surfactants. 

The  development  of  the  solutions  to  (2)-(8)  beginning  with  the  outer  solution  for  ^2  is  now 
straightforward.  The  matching  condition  for  the  second-order  pressure  gives  the  group  velocity 
Cf .  At  third  order,  where  we  terminate  the  analysis,  the  matching  condition  yields  an  evolution 
equation  for  A((,  r)  in  the  form  of  a  dissipative  NLS  equation  that  includes  the  effects  of  viscosity, 
surface  tension,  and  surfactant.  This  dissipative  NLS  equation  can  be  used  to  obtain  the  damping 
coefficient  and  to  examine  the  stability  of  weakly  nonlinear  waves. 
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3  Gravity  wave 


We  first  examine  the  case  when  the  capillary  effect  is  absent;  i.e.,  we  take  T  =  0  and  then, 
by  necessity,  k  =  0.  The  wave  amplitude  is  chosen  to  be  of  the  same  order  as  the  boundary-layer 
thickness,  so  that 

Aa  =  ^l/^  (34) 

where  A  =  0(1)  is  a  proportionality  constant.  The  gauge  functions  for  the  inner  expansions 
should  then  be 

Mc)=(— j.  (35) 

where  A  is  inserted  for  convenience.  We  have  verified  at  each  step  that  this  is  the  correci,  ^^equence 
for  the  first  few  terms.  The  first-order  inner  and  outer  solutions  are  identical  to  those  in  the 
previous  section  except  that  G  is  replaced  by  unity,  and  Qx  and  Qo  in  (31)  and  (32)  are  absent; 
thus  uj  is  independent  of  y* . 

The  second-order  term,  0(c),  in  the  potential  is  found  from  (13), 

<t>'2xx  +  4>2yy  —  —  2iA(Ee^  -h  c.c.  ,  (36) 

and  we  set 

=  iA^E^e'>  -  iAfEye*  +  c.c.  +  r).  (37) 

The  solution  to  the  homogeneous  equation,  fj(^,  y,  T)E^e^+Fi{^,  y,  T)Ee*+c.c.,  is  not  included. 
Matching  shows  that  Fj  =  0.  The  other  coefficient  Fi  remains  arbitrary  and  may  be  set  to  zero; 
taking  F^  ^  0  would  be  equivalent  to  changing  .<4  by  an  amount  aFi.  The  long-scale  velocity 
potential  should  satisfy 

=  0>  (38) 

which  can  be  obtained  by  extending  the  expansions  (12)-(14)  to  the  fourth  order.  (FVom  another 
viewpoint,  the  flow  at  large  distances  such  that  ^  =  0(1)  and  y  =  0(1)  is  described  by  (38).) 
The  boundary  condition  for  (^5  will  be  determined  at  third  order.  The  second-order  velocity 


1 

1 

1 

components  and  the  pressure  in  the  outer  region,  from  (17)-(19),  become 

1 

U2  =  A{E(e»  +  ye")  -  +  c.c.  +  2  |  Ap  e" 

(39) 

1 

Vi  =  iA^E^e* -iA(E{e« +  ye»)  +  c.c. 

(40) 

1 

P2  =  -A^E^e"  +  A(E{cge*  +  ye")  +  c.c.  +  2  |Ap  (e"  -  e^"). 

(41) 

1 

In  the  boundary  layer,  the  second-order  equations  are 

1 

1  . 

J«2v*  =  -«i. 

(42) 

■ 

»72i  =  uj  -  uJjji*  +  CgTiif  at  y*  =  0 

(43) 

1 

1  . 

jP2y*  —  “Wl* 

(44) 

1 

pj  =  »J2  at  y*  =  0 

(45) 

■ 

-  «2t  =  P2*  -  J»?l*P2y  +  P*(  - 

(46) 

1 

— «2y.  =  at  y  =  0. 

(47) 

■ 

After  matching  with  the  outer  solutions,  the  solutions  to  (42)-(47)  become 

■ 

uj  =  AA£y*  +  iA^E^  -  iA(E  +  c.c. 

(48) 

1 

^2  =  —A^E^  -f  ^A(E  -1-  c.c. 

4t 

(49) 

1 

p5  =  iXAEy*  -  A^E^  +  ^A(E  +  c.c. 

(50) 

m 

«;  =  A>/2(1  -  i)  exp  (^y*)  AE  -b  iXAEy*  -  A^E^  +  A(E  +  c.c.  +  2  |  A 

(51) 

1 

A  long-scale  free-surface  elevation  r^{^,r)  that  could  appear  in  (49)  as  a  result  of  integrating 

(43)  has  been  set  to  zero  as  a  consequence  cf  matching  the  pressure.  Also,  the  group  velocity  c. 

1 

has  been  replaced  by 

i 

The  third-order  outer  solutions,  0(«^/*),  are  obtained  first  by  solving 

1 

^3»»  -b  4>3py  —  12A^.E*e"  —  AAfE^{'^c*  -b  3ye*) 

1 

-AAJ(|e»  -b  ye")  -  A«^(e»  -b  2ye»)  -b  c.c. 

(52) 

1 

1 
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The  solution  to  the  homogeneous  equation  is  dropped  for  the  same  reason  as  cited  for  ^  in  (37). 
We  try  a  solution  for  (52)  of  the  form 

^3  =  +  AA^E\\ty  +  ye*)  -  AX|(|e»  +  ye»)  -  ^A^Ey^ey  +  c.c.  +  ^°(^ ,  y,  r).  (53) 

The  outer  solutions  for  the  vertical  component  of  the  velocity  and  the  hydrodynamic  pressure 
are  then  obtained  from  (18)-(19); 

U3  =  -|A3£;V  +  AA<£;»(|e*  +  ye«')-AA|(|e«'  +  yc») 

-\A((E(2  +  y)yc»  +  c.c.  +  (54) 

P3  =  -i|  +  iAA|(e»  -  e’»  -  ye**)  -  |A((r(l  +  y)ye* 

-ArEey  +  i  I A  I*  A£;(3e*  -  4e*»)  +  iAA^E^^e*  +  e*»  +  ye*»)  +  c.c.  +  (55) 


The  boundary-layer  equations  at  third  order  are 

“  “a*  “  (56) 

V3t  =  vl  +  ^Ti2(  -  Tiir  -  uj(»;2*  ■+■  r)i()  -  u^rjir  at  y*  =  0  (57) 

=  ~*’2I  +  2’'u  “  (56) 

pI  =  2AwJ,,  -b  ij3  at  y’  =  0.  (59) 


The  equations  for  U3  are  omitted  since  they  are  not  required  for  derivation  of  the  desired  aixq>li- 
tude  equation.  After  matching  with  (18),  equation  (56)  gives 

vl  =  -\s/2{\  -  i)exp  1^1’)  +  AiA*f?*y* 

-|-A*iA£:y**  -  2iAAf  Ey*  -  |a*E*  -b  fAA^E*  -  |aA|  +  c.c.  -b  4>\,.  (60) 

Integrating  (57)  then  yields  the  third-order  free-surface  elevation: 

>73  =  -2»A*E*  -b  3t  I A  I*  AE  -b  2A*AE  -b  2tAA<E*  +  ^A«E  -b  ArE  -b  c.c,  -b  >7|(^,  r),  (61) 
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where  the  condition 


A  I  4  la 

=  at  y  =  0  (62) 

has  been  imposed  to  suppress  secular  terms.  The  equation  for  pj  with  the  boundary  condition 
(59)  can  easily  be  solved  to  yield 

p;  =  ^X^AEy'^-XA^E^y’  +  ^XA(Ey‘-hAX^AE-^iA^E^  +  3i\A\^  AE 

2iAA(E^  +  ^A((E  +  ArE  +  c.e.  -  2A  |  A  I*  y*  +  q§(^,  r).  (63) 

The  matching  condition  for  the  third-order  pressure  (63)  with  the  outer  solution  (19)  finally  gives 

at  »  =  0.  (64) 

for  terms  independent  of  E,  and 

iAr  +  2iX^A-lA((  =  2\A\‘^  A,  (65) 

Q 

for  terms  linear  in  E. 

Except  for  an  additional  term  2tA^i4,  the  amplitude  equation  (65)  is  identical  to  the  NLS 
equation  that  governs  the  amplitude  modulation  of  gravity  waves  in  inviscid  flow.  It  is  thus 
appropriate  to  call  the  amplitude  equation  a  dissipative  NLS  equation.  ^FVom  (65),  the  integral 
of  |v4 1*  over  ^  is  easily  found  to  decay  as  exp(-4A*T). 

For  infinite  wavetrains  independent  of  the  solution  Ao(r)  of  (65)  that  satisfies  i4o(0)  =  ao 
is 

^o(r)  =  aoexp  +  ^  l®ol*  (e"**’’’  -  l)j  .  (66) 

where  uq  is  a  complex  constant  with  |ao  |=  1/2,  as  explained  earlier.  Here,  the  expression  |ao 
is  kept  to  show  the  amplitude  dependence  explicitly.  In  the  inviscid  limit,  A  — >  0,  the  Stokes 
wave  (to  third  order)  is  recovered: 

Ao{t)  =  no  exp  (-2i  |  oo  |*  r)  .  (67) 
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The  decaying  factor  in  (66),  combined  with  (9)  and  (34),  becomes  exp(-2d).  Thus,  the  damping 
coefficient  Dg  for  weakly  nonlinear  gravity  waves  is 

D,  =  2€,  (68) 

which  is  identical  to  the  linear  result  of  Stokes  (1845).  In  retrospect,  we  could  instead  have 
shown  that  the  damping  coefficient  is  unaffected  by  nonlinearity  at  this  order  using  the  dissipation 
function  and  its  expansion  for  small  amplitude.  Viscosity  also  causes  a  phase  shift,  as  can  be  seen 
from  (66).  Although  viscous  dissipation  causes  exponential  decay  in  the  amplitude  of  nonlinear 
waves,  it  is  not  clear  whether  it  can  suppress  a  modulational  instability  of  the  Benjamin-Feir 
type  and  the  subsequent  recurrence  phenomenon.  We  will  examine  this  in  §5. 

4  Capillary-gravity  waves 

We  now  consider  capillary-gravity  waves  with  a  surfactant  on  the  free  surface.  For  the  gravity 

0 

wave  in  the  previous  section,  the  amplitude  of  the  wave  is  chosen  to  be  of  the  same  order  as 
the  boundary-layer  thickness  so  that  nonlinearity  and  dissipation  both  appear  in  the  amplitude 
equation.  When  a  surfactant  is  present,  however,  lower-order  dissipation  is  expected;  thus,  the 
amplitude  is  chosen  to  be  of  lower  order  than  the  boundary-layer  thickness: 

Ao  =  (69) 

where  A  is  again  a  proportionality  constant.  The  gauge  functions  for  the  inner  expansions  now 
become 


The  first-order  solutions,  are  as  given  by  (17)-(19)  for  the  outer  region  of  irrotational 

flow  and  by  (28)-(32)  for  the  boundary  layer.  Since  the  nondimensional  surface  dilational  modulus 
K  and  the  Weber  number  T  are  nonzero,  the  boundary-layer  solutions  for  the  h<»izontal  velocity 
component  (31)  and  displacement  (32)  are  affected  by  the  surfactant  at  the  leading  order,  and 
the  linear  dispersion  relationship  (33)  is  that  for  capillary-gravity  waves. 
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The  equation  for  the  term  4>i  in  the  velocity  potential  is  identical  to  (36),  but  the 

appropriate  solution  now  includes  a  solution  of  the  homogeneous  equation: 

(t>2  =  -  iA(Eye*  +  c.c.  +  ^§(4,  y,  r),  (71) 

where  F(^,  r)  can  be  determined  through  a  matching  condition.  Again,  the  long-scale  potential 
^2  should  satisfy  the  Laplace  equation  (38),  with  the  boundary  condition  to  be  determined 
through  matching  of  the  third-order  solutions.  The  velocity  components  and  the  pressure  in  the 
outer  region,  from  (17)-(19),  then  become 

«2  =  2iFE^e^^ +A(E(e» +  ye*)-A^E^e'> +  c.c.  +  2\A\^  (72) 

V2  =  2FE^e^^ +  iA^E^e» -iA(E(e» +  ye^)  +  c.c.  (73) 

P2  =  2iFE^e^^ -A^E^e* +  A(Eie,e> +  ye»)  +  c.c.A2\A\^  {e* -e^).  (74) 

The  boundary-layer  equations  for  the  0(€*^*)  terms  differ  from  those  for  the  gravity  wave  in 
the  previous  section  because  the  amplitude  considered  here  is  of  lower  order  than  the  boundary- 
layer  thickness.  Substituting  the  inner  expansions  into  (2)-(8)  and  using  the  multiple  scales  (9), 
the  inner  variable  (20),  and  the  amplitude  (69),  we  obtain 


(75) 

TUt  =  vl  -  ulrii,  +  CfftK  at  y*  =  0  (76) 

rf,.  =  0  (77) 

Pi  —  Offi  —  T{ij2xt  +  at  y*  =  y  (78) 

”  ’*5*  ~  rir  +  Pl(  ~  +  “I “I*  “  (79) 

«5f  =  ''{(Cs  +  »7i>7i<r),,  +  2Ci«<}  at  y*  =  0.  (80) 


The  third  ^proximation  to  the  y*  momentum  equation  is  also  needed,  to  elinninate  from 
(79): 


(81) 
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The  vertical  velocity  component  V2  is  obtained  by  integrating  (75)  and  imposing  the  matching 
condition; 

=  -<3i  I A  p  +2FE^  +  (1  -  Qi)A^E^  -  iA(E  +  c.c.  (82) 

Then,  from  (76),  the  second-order  free-surface  elevation  becomes 


Ti2  =  iFE^  -  A^E^  +  (1  -  c,)A(E  -I-  c.c.,  (83) 

where  the  long-sc2de  quantity  f;3(^,y,r)  is  omitted  in  anticipation  of  the  matching  condition  for 
the  second-order  pressure.  After  (83)  is  substituted  into  the  boundary  condition  for  pj  in  (78), 
the  inner  solution  for  the  second-order  pressure  is  obtained: 


pj  =  (l-t-  ZT)iiF  -  A»)r*  H-  (1  -  c,  -f  2T)A(E  +  c.c.  (84) 

The  solution  (84)  should  match  the  outer  solution  (19),  with  the  terms  0(f*/^)  replaced  by  (74). 
Therefore,  the  slowly  varying  rmiplitude  F  and  the  group  velocity  Cg  are  obtained  as 

(85) 

(86) 


i+2r 
2  ’ 


Ca  — 


respectively.  The  group  velocity  (86)  agrees  with  the  relationship  for  a  linear  capillary-gravity 
wave  in  inviscid  flow.  It  is  noteworthy  that  F  blows  up  as  T  — »  1/3;  thus,  a  second-harmonic 
resonance  persists  despite  the  effects  of  viscosity  and  surfactants.  A  complete  discussion  of  this 
resonant  behavior  was  given  by  McGoldrick  (1970).  A  modified  asymptotic  flow  description  is 
needed  for  small  values  of  |T  —  1/3|  and  will  be  discussed  later  in  this  section. 

Substituting  the  known  solutions  into  the  right-hand  side  of  the  boundary-layer  equation  for 
U]  in  (79)  gives 


tijy.y.  —  Uji  = 


iQl  _  2Qi  -  A^E^  +  (•  -  CgQy)A(E  +  c.c. 


(87) 


The  boundary  condition  for  uj  in  (79)  can  be  ^>plied  more  conveniently  after  it  is  differentiated 
with  respect  to  time  and  the  kinematic  boundary  condition  (8)  is  introduced: 


~  ««j*«  =  -4»«(Qo  +  2tQo  -  3)A^E’  +  k{2  -  Cf){iQo  -  1)A(E  +  c.c.  at  y*  =  0.  (88) 
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Then,  from  (87)  and  (88),  the  second-order  horizontal  velocity  component  is 


-  \q\  +  2i(5ij  +  |?2exp  ^ 


l-f-ST 

i-3r 


+  c.c  +  2\A\\ 


where 


=  -2K-\-i  [(^  +  «)  <?2  +  {V5+  ‘•(8'C  -  V2)}Qo  +  4«|-||] 


amd 


^  [i2v/2(2  -  c,)  -K1  -  s)cj  Qo  -  2v^(l  -  c,) 

2(\/5k-1-») 

In  the  third  approximation,  the  velocity  potential  ^3  should  satisfy  (14),  or 

-AfiEie^  -I-  2ye»)  -h  AA(E^  -  3ye»  + 

-AAl  ^^^~e»  +  ye»)d-c.c. 


(89) 


(90) 


The  appropriate  solution  of  (90)  is 


«3  = 

+AA,E’  +  i"*) 

+c.c.-»-^3(^,y,T),  (91) 

where  the  slowly  varying  amplitudes  f3((,  r)  and  F3((,  r)  are  to  be  determined  from  the  matching 
condition.  As  in  the  derivation  of  ^3,  a  first-harmonic  term  in  the  solution  to  the  homogeneous 
equation  is  drc^ped.  The  vertical  velocity  component  V3  is  then  obtained  from  (18)  as 


V3 


=  +  2F.r>.«.  -  A’E’  +  M|>  AE^, 


..  .  r,a /'5-2T  .  ,  6r  3,  12r  -.X  /'3-2T. 

+AA(E^  — - — e*  +  ye»  +  ^  _  gy*  *  +  )  ~  — 2 — ®  j 

~A„F  (ye»  +  ^e*)  c.c.  (92) 
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Replacing  the  terms  in  (19)  with  the  known  solutions  yields 


P3  =  +  i2F2E^e^*  -iA^E^ 

+iAA,E^ 


.  iA^E^  ^  _  1-^“^  _  _«2_,3, 

2(1 -3T)  1-3T  1-3T 

i-er-er^ 

1  -  3r  1  -  37’  ^ 


+.■  \AfAE  [  -  ^ew‘"’  -  tS'”]  -  (P  + 

+ii4j4|(e*  —  e^*'  —  ye^)  —  ArEe^  +  e.c.  +  (iju  —  4>iy)<^2y  +  (93) 

The  third-order  boundary-layer  solutions  are  again  obtained  from  the  continuity  equation  (4) 
for  uj,  the  free-surface  boundary  condition  (7),  and  the  vertical  momentum  equation  (3)  for  pj. 
The  equation  for  the  vertical  velocity  component  then  becomes 


*^3y*  —  ~^***l*  +  ’7l***2y»  +  ’ISrUly*) 


which  gives 


=  [(^Qi+>^V)AE  +  c.c]-Y^iiA^E^  +  c.c){QiAE  +  c.c)-{AE  +  c.c)- 


-  ig?  +  2iQ,)A^E^  +  {5,exp  (^v’)  -  +  c.c 

1  +  7T  1  1  —  YT 

+  [3{f3  -  2(i-3r)^^^^  +  H  ^  +  c.c)(giyif;  -i-  c.c) 

+  ^AfAE-i^AAt+  ^AAtE^  A  c.c]  (95) 

after  the  matching  condition  is  applied.  The  third-order  free-surface  elevation  73  is  obtained  by 


integrating 


»J3«  =  t'3+Cr 'Tat -*71 T -«!(%,-»- »7i()-u;?7i,  at  y*  =  0. 


The  result  after  appropriate  substitutions  is 


+  1  ‘  ^4  ~  +  c.c.  +  773  (97) 


where  the  condition 


^5,  =  2(1  -  T)  |<  -iiQoAAl  -  QIA^A() 


at  y  =  0 
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has  been  imposed  to  suppress  secular  terms.  As  mentioned  earlier,  (98)  is  the  boundary  condition 
for  ^2!  thus,  ^2  ^  uniquely  determined  once  A  is  obtained  from  (104)  below.  It  is  noteworthy 
that  additional  nonlinear  forcing  terms  are  present  in  (98)  because  of  the  surfactant. 

The  differential  equation  for  P3  is  (81),  and  the  associated  boundary  condition  is 

P3  =  +  Grfi  —  T(rfy„  +  (99) 


Solving  (81)  subject  to  (99)  results  in 

...... 

.(1  +  3r)(4-  i3r+6T^)  +  ler  i  +  3r 


P3  =  « 

+ 

+i 


2(1  -  3T) 

.l-8T  +  4r2  ,  „  ,  ^  . 

Ar .E  +  t 


+  — 


AA^E^  +  A*(iy  +  I^Q^)AE 


.6  +  3T  +  9T2 


4  '  ■  2(1 -3T) 

Applying  the  matching  condition  for  the  pressure,  we  determine  and  Fz: 


(100) 


37(5  4- 2ir)  a 
4(1  -  3T)(1  -  4T)  ’ 


(101) 


'37(7-157  +  67^) 
2(1-37)2 


■  1  +  37 
*2(1  -  37) 


AA^. 


(102) 


We  also  obtain 


»73 


1  +  27 
2(1  -  7) 


Ai 


(103) 


and  the  amplitude  equation  for  the  capillary-gravity  wave  with  viscosity  and  surfactant: 


»At  + 


iA2 
2^  ' 


k[k  +  i(\/2  -  k)]  ^ 


i(4T>-8r+lM(f 


972  -  157  +  8 
4(1  -  37) 


|A|>  A. 


(104) 


Two  more  nonuniformities  are  observed,  in  addition  to  the  one  for  the  second-harmonic  res¬ 
onance  (7  =  1/3).  In  (101),  the  amplitude  F3  becomes  unbounded  as  7  approaches  1/4;  thus, 
a  third-harmonic  resonance  occurs.  As  7  approaches  unity  (the  value  for  pure  capillary  waves), 
^3  in  (103)  is  singular,  and  a  rescaling  of  the  Itmg-scale  free-surface  elevations  is  required,  as 
explained  below. 
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The  dissipation  term  proportional  to  A  in  (104)  is  attributed  to  the  surfactant,  whereas  for  a 
clean  surface  the  dissipation  is  of  higher  order,  as  shown  in  the  previous  section  where  different 
scaling  was  required.  Therefore,  if  either  A  or  k  is  zero,  the  dissipation  term  in  (104)  is  absent, 
and  we  recover  the  amplitude  equation  for  capillary-gravity  waves  in  inviscid  flow.  In  these  cases, 
the  resulting  NLS  equation  could  also  be  obtained  more  easily  in  a  Cartesian  coordinate  system 
by  incorporating  the  vertical  modulation  y  in  the  derivation  of  Djordjevic  ic  Redekopp  (1977) 
for  finite  depth.  In  the  dissipative  NLS  equation  (104),  the  linear  dissipation  term  has  a  complex 
coefficient;  the  imaginary  part  is  related  to  the  decay  rate,  and  the  real  part  corresponds  to  a 
frequency  change  due  to  the  surfactant. 

A  solution  of  (104)  for  infinite  wavetrains  which  depends  only  on  r  is 

2  f9T^-15T  +  8\ 


Ao  =  ao  exp 

f  I 


A2  /c2-(-i(V5k-/c2)  ‘ 

2>/2  /c2-v^x-j-l 

•  exp 

*-|ao|= 


V  1-3T 


(105) 


2v^/c2A2 

Since  the  slow  time  is  now  defined  as  r  =  due  to  (69),  the  damping  coefficient  D,  for 

an  infinite  wavetrain  (105)  becomes 

D,  =  ^  (106) 

'  4(k»- V2IC-H)  ^  ’ 

which  has  a  maximum  at  k  =  \/2.  At  this  maximum,  the  frequency  shift  due  to  the  surfactant 

disappears.  The  damping  coefficient  (106)  agrees  with  the  results  of  Levich  (1941,  1972)  amd 

Lucassen-Reynders  k  Lucassen  (1969)  for  linear  capillary  waves,  as  can  be  seen  by  expanding 

their  solutions  for  large  Reynolds  number.  It  is  clear  from  the  result  that  the  damping  effect  is 

greatly  enhanced  by  the  surfactant. 

In  the  derivation  of  these  expansions,  we  have  noted  that  singular  behavior  appears  for  certain 
values  of  the  Weber  number.  The  case  for  T  =  1/3  corresponds  to  second-harmonic  resonance  and 
is  related  to  Wilton’s  ripples.  The  nonuniformity  can  be  removed  by  considering  a  superposition 
of  the  fundamental  wave  and  its  second  harmonic  in  the  first  approximation  and  by  adding  an 
intermediate  slow  time  scale,  as  in  McGoldrick  (1970)  for  inviscid  flow.  We  start  by  considering 
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a  Weber  number  near  the  value  for  the  resonance; 


where  T  =  0(1)  is  a  constant.  An  appropriate  solution  for  (12)  is  now 

<t>i  =  Ai(4,f)£e»  +  +  c.c., 


(107) 


(108) 


instead  of  (15),  where 


i=at  =  r/a 


(109) 


is  the  new  time  variable.  The  subsequent  analysis  proceeds  as  before  except  for  the  added 
complexity  due  to  the  presence  of  second-harmonic  terms  in  the  iirst-order  solutions.  For  example, 
an  equation  analogous  to  (36)  is 


+  =  -24tAifr‘e*»-ii4,Aj£®(10e*»-|-8e»)-3«A?EV 

-6»AJ -  2iAi(Ee»  -  AiAifEe^*  +  c.c.,  (110) 


and  an  appropriate  solution  of  the  homogeneous  equation  should  include  third-  and  fourth- 
harmonic  terms.  A  matching  condition  for  the  seccmd-order  pressure  finally  gives 


Au  =  -AIA2  (111) 

TT  1  1 

A2i  +  iYA2  +  -^A2(  =  -Al  (112) 

in  addition  to  the  group  velocity  c,  =  5/6,  the  ccM’rect  value  when  T  =  1/3. 

The  evolution  equations  (111)-(112)  are  equivalent  to  the  second-harmonic  resonance  equa¬ 
tions  given  by  McGoldrick  (1970),  except  that  a  term  proportional  to  7*  is  added  to  account  for 
values  of  the  Weber  number  close  to  1/3.  The  coefficients  for  the  terms  with  spatial  variation 
are  different  because  the  slow  variable  (  is  constant  at  points  moving  with  the  group  velocity; 
adding  Cgd/d^  to  d/^ leads  to  the  same  terms  as  in  McGoldrick  (1970).  For  large  |T|,  equations 
(111)-(112)  are  consistent  with  previous  results.  As  ITj-^  oo,  (112)  reduces  to 


which,  with  (107),  is  found  to  agree  with  (85).  Using  (113)  to  remove  A2  in  (111)  gives 

Mil*  ^1-  (114) 

This  equation  can  also  be  obtained  from  (104)  by  using  (107),  setting  r  =  at,  and  noting  that 
the  dissipation  and  the  dispersion  terms  become  small  as  T  1/3.  Thus,  the  results  for  T  ^  1/3 
and  for  T  close  to  1/3  match  asymptotically  as  T  —►1/3  and  [T-  1/3|  /o  — ►  00. 

A  set  of  solutions  of  (111)-(112)  for  infinite  wavetrains  without  (-dependence  can  be  obtained 

as 

Ai  =  6oe~‘**  (115) 

’  (I) 

where 

6  =  \{it-  (sgnf)v^9f’-H6|6o|»)  ,  ( 117) 

and  60  is  a  complex  constant.  The  sign  in  (117)  has  been  chosen  so  that  0  -*  0  when  ITf-*  00. 
Then,  also  1  Aj  |— •  0,  and  hence  60  becomes  identical  to  oo  as  j  T  |—  00.  When  f  is  zero,  the 
wavetrains  obtained  by  McGoldrick  (1970)  are  recovered  from  (115)-(117). 

A  nonuniformity  related  to  the  third-harmonic  resonance  is  observed  when  T  =  1/4.  We 
can  expect  more  nonuniformities  for  T  =  l/(n 1)  (n  =  3,4,  -).  corresponding  to  a  n‘*- 
harmonic  resonance.  Uniform  expansions  near  these  resonance  values  of  the  Weber  number  can 
be  obtained  by  modifying  those  introduced  for  T  —  1/3  above.  The  amplitude  equation  (104)  still 
can  describe  the  evolution  of  the  slowly  varying  amplitude  of  the  first-order  solutions,  because 
higher-harm<»ic  corrections  need  not  be  imposed  on  the  first-order  solutions  (Joo  et  al.,  1989). 

Still  another  non<'nifonnity  occurs  for  values  of  T  close  to  unity,  when  the  surface-tension 
force  is  large  in  comparison  with  the  gravitational  force.  In  this  case,  the  long-wave  solutions  for 
the  velocity  potential  and  the  surface  elevation  must  be  modified.  Since  the  orders  of  magnitude 
change,  it  is  now  convenient  to  omit  subscripts  and  to  denote  the  largest  “long-wave”  terms 
by  and  where  4^  and  1;®  are  not  necessarily  0(1).  By  repeating  the  previous 
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derivations  when  1  —  T  <  1,  we  find  that  the  equations  analogous  to  (98)  and  (103)  become 

-  ~  i(QoAAt  -  QlA‘Ai)  (118) 

(l_r),,°  =  i-!^ci/V{°.  (119) 

When  1  —  r  =  0(1),  then  =  0(1)  but  i]°  =  0(f*/^);  t)°  and  <t>°  can  be  replaced  by  and 
^§,  respectively,  and  (98)  and  (103)  are  recovered.  When  <;  1  —  T  C  1,  »7°  =  0{e^/^/(l  — 
r)}  while  remains  0(1).  In  the  distinguished  limit  corresponding  to  1  —  T  =  0(e^/'*),  all 
terms  in  (118)  and  (119)  are  retained.  Finally,  for  1  —  T  <C  rj°  —  0(1)  while  = 

0({(1  -  T)/(}l*)).  Thus,  for  pure  capillary  disturbances  in  the  presence  of  a  surfactant,  there 
is  a  long-wave  component  with  a  surface  elevation  whereas  the  corresponding  long¬ 

wave  potential  is  of  higher  order;  r;°  is  determined  by  (118).  For  a  clean  surface,  however,  the 
nonuniformity  at  T  =  1  disappears  because  Oo  =  0  and  is  proportional  to  1  —  T. 

5  Linear  stability  analysis 

The  effects  of  viscosity  and  surfactant  on  the  stability  of  infinite  wavetrains  to  sideband 
disturbances  of  the  Benjamin-Feir  type  can  be  examined  by  considering  the  dissipative  NLS 
equations  (65)  and  (104).  We  can  obtain  the  evolution  of  disturbances  to  the  solutions  (66)  and 
(105)  for  small  time,  while  the  magnitude  of  the  disturbances  stays  relatively  small,  by  extending 
the  linear  analysis  of  Stuart  k  DiPrima  (1978);  this  extension  incorporates  the  unsteadiness  of 
the  basic  state. 

We  first  consider  a  perturbation  to  the  wavetrain  (66)  for  gravity  waves  and  write 

A  =  Ao(r)ll-hB({,T)].  (120) 

Substituting  (120)  into  (65)  gives 

iB,-iB„=21o<,|’e-«*'(B-*-B*)  (121) 

after  the  quadratic  ternw  in  B  are  neglected.  For  a  pair  sideband  modes,  we  seek  a  solution 
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of  (121)  in  the  form 


(122) 


where  Bi  and  B2  are  complex  coefficients,  and  /  is  a  real  wavenumber. 

Although  the  disturbance  can  grow  for  small  values  of  the  slow  time  variable  r,  it  will  even¬ 
tually  be  damped  out  due  to  viscous  dissipation.  A  norm  for  one  spatial  period, 


N{t)  =  I  |A|  -  |Ao| 


(123) 


can  be  used  as  a  measure  of  the  growth  or  decay  of  the  disturbance.  While  the  magnitude  of  the 
disturbance  B  stays  small,  the  norm  N  is  equivalent  to 


N(r)  =  -^  (|Bi  |2  +  |B,  I"  +B^B2  +  B\Blf'^  ,  (124) 

which  is  seen  by  substituting  (120)  and  (122)  into  (123). 

We  now  substitute  (122)  into  (121)  and  obtain  a  pair  of  linear  homogeneous  equations  for  Bi 
and  B^- 


5l--jBi  +  2i|ao|'e-‘^*"(Bi  +  BJ)  =  0 

fli-yfl2-b2i|ao|’e-‘^’"(Sf-|-B2)  =  0, 


(125) 

(126) 


where  a  prime  denotes  differentiation  with  respect  to  r.  Equations  (125)-(126)  and  the  transfor¬ 
mation 

fli  =  e-’**^B(r)  (127) 


yield 

=  (128) 

where 

^ (32  I  ao  P  e--***"  -  f»)  +  I +  ( 129) 

In  the  inviscid  limit  (A  -*  0),  |  Ao  |  is  constant,  and  (129)  becomes 

«  =  ^(32  loop-/*)  (130) 
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and  B  =  Bi-  The  stability  of  the  flow  is  then  determined  by  the  sign  o{$:  0  >  0  indicates  insta¬ 
bility,  whereas  d  <  0  gives  neutral  stability,  in  agreement  with  the  inviscid  instability  condition 
of  Benjamin  k.  Feir  (1967): 

0</^<32|ao|*.  (131) 


The  details  of  the  instability  for  inviscid  flow  were  given  by  Stuart  k  DiPrima  (1978). 

When  viscosity  is  present,  the  solution  for  B  is  obtained  from  (128).  A  change  of  variables 


_  |aoM  ^-aA»T 

2V2A 


(132) 


leads  to  the  Bessel  equation 


where 


The  solution  for  B  is 


^d}B  ,  dB 


a  o  ao  ,  t  5.  = 

+  +a)B  =  0, 


a  =  1  -1- » 


16A2‘ 


B  =  Cl  A’o(s)  +  C2/a(s), 


(133) 


(134; 


(135) 


where  Ka  and  are  modifled  Bessel  functions  of  complex  order  a,  and  ci  and  C2  are  constants 
to  be  determined  by  initial  conditions  at  s  =|ao|  //(2'\/2A*).  For  large  t,  and  thus  small  s,  the 
disturbances  are  expected  to  decay  as  a  result  of  viscous  damping.  As  5  — »  0,  the  Bessel  functions 
in  (135)  have  the  form 


~  (constant)s~®,  /o(s)  ~  (constant)s“.  (136) 

The  solution  for  B,  therefore,  is 

B  ~  (con8tant)e“*’^/*e’*’'  (137) 

as  r  — ►  00.  However,  since  Bi  =  e“^**’'5  and  |  Aq  |=  ooc"***’’,  the  amplitude  of  the  disturbance 
as  T  — ►  00  is 

|Aofl,|=0(e-2**")  (138) 
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and  approaches  zero  as  anticipated. 

The  result  (138)  obviously  does  not  allow  recovery  of  the  condition  (131)  obtained  for  zero 
viscosity.  If  A  is  small,  both  the  order  and  the  argument  of  the  Bessel  functions  are  large  in 
magnitude.  To  achieve  a  single  asymptotic  representation  that  contains  (131)  and  also  predicts 
decay  for  large  time,  we  take  the  limit  as  A  — ►  0  and  t  — ►  oo  with  A^r  fixed.  Then, 


Us)  = 


2  \  -1/-* 


exp  +  _  +  ^=_ 


The  previous  results  are  now  recovered  for  large  and  small  values  of  A^r.  If  A^r  — ►  oo,  the  largest 
term  arises  from  Ka  and  is  proportional  to  8“®,  as  in  (136);  thus,  B  is  again  found  in  the  form 


(137).  If  A^r  —  0, 


i.  =  |ooM 

a  2v^A2a 


(l>2A2r+.), 


and  B,  after  some  manipulation,  becomes 


B  (constant)  exp 

+ (constant)  exp  y^l  -  — ■  .  (142) 

This  is  exactly  the  solution  of  ( 128)  when  A  =  0. 

The  linear  stability  analysis  for  capillary-gravity  waves  represented  by  (105)  proceeds  in  the 
same  way  as  for  the  gravity  wave.  The  evolution  equation  for  the  disturbance,  equivalent  to 


(121),  becomes 


The  transformation  and  definition  for  9  corresponding  to  (127)  and  (130)  are 

=  e-®**^B(r)  (144) 

6  =  Ti/^  (2  |ao  I*  -  Ti/*)  +  i2£>TiA®/=  •¥  D^X*,  (145) 
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respectively,  where 


D  = 


2^  ic*- v/2K  +  r 


and 


47’2_8r+l  .  ^  9T2-15r+8 

Ti  = - r -  and  Tj  = 


8 


4(1  -  3T) 


For  inviscid  flow,  (145)  reduces  to 

e  =  Ti/2  (2  |ao|^  r,  -  Ti/*) 


(146) 

(147) 

(148) 


When  TiTj  <  0,  or 


(149) 


6  is  always  negative,  and  the  side-band  resonance  disappears  regardless  of  the  disturbance 
wavenumber.  The  corresponding  stability  analysis  for  waves  in  water  of  finite  depth  has  been  re¬ 
ported  by  Pjordjevic  b.  Redekopp  (1977).  The  stability  band  (149)  corresponds  to  the  deep-water 
limit  of  their  stable  region,  which  lies  between  the  minimum  phase  speed  and  the  second-harmonic 
resonance.  The  other  stable  regions  seem  to  disappear  as  the  depth  approaches  infinity.  Although 
only  neutral  stability  is  predicted  for  (149),  computations  of  the  large-time  modulation  for  this 
region,  in  the  following  section,  show  strikingly  different  behavior  from  that  in  the  other  neu¬ 
trally  stable  regions.  When  0  <  T  <  1  -  ■\/3/2,  the  lower  instability  boundary  for  /  stays  at  sero, 
whereas  the  upper  bound  increases  with  T  until  it  blows  up  at  T  =  1  —  y/i/2.  When  the  Weber 
number  exceeds  the  value  for  the  second-harmonic  resonance,  the  upper  bound  decreases  with 
the  Weber  number.  As  T  — »  1,  the  upper  bound  approaches  /  =  1/v^,  which  is  smaller  than  the 
value  /  =  2‘)/2  for  gravity  waves. 

For  viscous  flow,  the  solution  for  B  with  (145)  should  be  considered.  The  change  of  variables 
analogous  to  (132)  is 


V2rira|ao|  f 


(150) 


The  resulting  equation  is  again  the  Bessel  equation  (135),  but  with 

.TiP 


o  =  1  -I-  I 


DA* 


(151) 
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The  subsequent  results  are  then  similar  to  those  for  the  gravity  wave  except  that  no  initial  growth 
is  possible  when  the  Weber  number  satisfies  (149),  as  anticipated.  A  large-time  modulation  for 
this  range  of  Weber  numbers  is  examined  in  the  following  computations. 

6  Numerical  solution  of  dissipative  NLS  equations 

The  linear  analysis  in  the  previous  section  is  valid  only  while  the  magnitude  of  the  disturbance 
remains  small.  In  this  section,  we  obtain  numerical  solutions  of  (65)  and  (104)  with  given  initial 
conditions  and  examine  the  evolution  of  a  disturbance  of  the  type  in  (122)  for  finite  time. 

Numerical  solutions  of  the  NLS  equation  have  been  obtained  by  Yuen  k  Ferguson  (1978) 
and  Weideman  k  Herbst  (1987)  and  many  others,  with  an  emphasis  on  modulational  instability 
and  recurrence  of  the  Fermi-Pasta-Ulam  type  (Fermi  et  rJ.,  1974).  Yuen  k  Ferguson  (1978) 
examined  the  relationship  between  Benjamin-Feir  instability  and  recurrence  and  showed  two 
types  of  recurrence,  “simple”  and  “complex.”  They  explained  that  for  modulations  with  the 
sideband  perturbation  wavenumber  /  in  the  range  8|ao  p  <  /’  <  32|ao  P  the  recurrence  is  simple 
because  all  higher  harmonics  of  the  prescribed  modulation  are  stable.  For  modulations  with 
0  <  /^  <  8|ao  p,  the  recurrence  is  complex  (quasi-periodic)  because  at  least  one  higher  harmonic 
of  the  prescribed  mode  lies  in  the  unstable  region.  Weideman  k  Herbst  ( 1987)  obtained  similar 
solutions  using  finite-difference,  spectral,  and  pseudospectral  methods,  and  they  compared  the 
effectiveness  of  these  methods.  We  use  a  pseudospectral  method  to  solve  the  dissipative  NLS 
equations  for  capillary-gravity  waves.  This  method  is  appropriate  for  our  purpose,  once  an 
aliasing  error  is  suppressed  by  introducing  a  sufficient  number  of  degrees  of  freedom  (number  of 
collocation  points). 

The  initial  condition  (120),  with  B  replaced  by  (122),  is  specialized  to  the  simple  case 

A  =  oo(l -t-4cos/(),  (152) 

where  the  substitution 

Bi(0)sB,(0)  =  ^6  (153) 
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has  been  made.  Modulational  behavior  for  different  initial  sideband  amplitudes,  Bi  and  Sj,  will 
be  discussed  later  in  this  section. 

A  periodic  boundary  condition 

A(^r)  =  AU  +  y,r),  (154) 

is  imposed  using  a  Fourier-coUocation  method  with  the  computational  domain  — x//  <  ^  <  x//. 

The  slowly  varying  amplitude  is  then  expressed  as 

N 

A=  5;  A„(r)e"‘<'/',  (155) 

n=-N 

where  we  choose  IV  as  a  positive  integer  power  of  2  and  a  fast  Fourier  transformation  is  used 
to  compute  the  Fourier  coefficients  An-  As  in  the  analyses,  the  amplitude  of  ao  is  and  6 
is  small.  In  particular,  we  set  oq  =  Ao(0)  =  0.5  and  6  =  2Ai(0)  =  2A_i(0)  =  0.1  for  most 
of  the  computations,  unless  otherwise  specified.  A  fourth-order  Hamming  modified  predictor- 
corrector  method  is  used  for  time  marching.  The  fourth-order  Runge-Kutta  method  is  used  for 
automatic  a4iustment  of  the  initial  tinoe  increment  and  for  computation  of  starting  values.  For 
most  computations,  the  maximum  time  step  is  0.01,  and  the  number  of  collocation  points  is  32 
for  one  initial  period  in 

Fig.  1  illustrates  the  evolution  of  the  Stokes  wave  (67)  for  zero  viscosity  when  subject  to 
the  modulation  described  by  (152).  Evdution  of  the  envelope  amplitude  |  A((,  r)  |  is  shown  in 
Fig.  la,  c,  and  e,  whereas  the  magnitude  of  each  Fourier  coefficient  |  An(r)  |  (n  =  0, 1,2,  ■  •  ■)  is 
plotted  in  Fig.  lb,  d,  and  f.  Since  |  An  |  is  symmetric  about  n  =  0  due  to  (152),  the  evolution 
f<»  negative  n  is  deleted  for  clarity.  The  perturbation  wavenumber  /  =  1  corresponds  to  complex 
recurrence  (Fig.  la,  b);  I  =  2  corresponds  to  the  maximum  initial  growth  rate  and  simple 
recurrence  (Fig.  lc,d).  These  figures  agree  with  the  results  of  Yuen  ic  Ferguson  (1978)  and 
Weideman  Sc  Herbst  (1987),  who  provide  detailed  explanations.  For  the  simple  recurrence,  the 
fundamental  (n  =  0)  and  sideband  (n  =  ±1)  inodes  are  periodic,  as  can  be  also  seen  in  Fig.  8. 
The  higher  inodes  (not  shown  in  Fig.  8),  excited  due  to  nonlinear  interaction,  are  not  exactly 
periodic,  but  they  are  not  dominant  at  any  stage,  making  the  modulation  almost  periodic  in  time. 
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When  /  lies  outside  of  the  instability  range  (131),  a  very  nearly  periodic  oscillation  is  observed 
(Fig.  le),  which  is  in  good  agreement  with  the  neutral  stability  predicted  by  the  linear  an^dysis. 
The  corresponding  Fourier-space  solution  (Fig.  If)  shows  that  none  of  the  higher  harmonics  is 
excited  and  that  the  sidebands  never  become  dominant. 

In  Fig.  2,  the  same  cases  as  in  Fig.  1  are  described  using  the  norm  N  defined  in  (123). 
The  behavior  for  small  r  is  in  precise  agreement  with  the  linear  analysis.  The  growth  rate  is  a 
maximum  when  /  =  2  and  decreases  as  /  increases,  resulting  in  larger  recurrence  periods,  until 
it  becomes  zero  for  /  =  2y/2.  When  /  becomes  larger  than  2%/2,  an  oscillation  due  to  neutral 
stability  is  observed,  as  predicted  by  the  linear  analysis.  When  /  decreases  from  2,  the  initial 
growth  rate  again  decreases,  but  complex  recurrence  is  observed  for  /  <  v^. 

Fig.  3  shows  the  evolution  of  the  slowly  varying  envelope  A  for  gravity  waves  with  small 
viscosity.  The  initial  condition  and  the  perturbation  wavenumber  /  are  identical  to  those  for  the 
cases  in  Fig.  1.  Here,  the  dissipative  NLS  equation  (65)  is  solved  with  A  =  0.125.  In  Fig.  3a 
and  3b,  the  values  of  /  lie  in  the  instability  range,  and  the  initial  behavior  shows  a  modulational 
instability  despite  viscous  dissipation,  which  is  in  agreement  with  the  analysis  in  the  previous 
section  for  small  A.  Compared  to  the  corresponding  cases  in  Fig.  1,  the  spikes  are  attenuated 
at  larger  times  for  both  the  complex  and  the  simple  recurrence.  In  Fig.  3c,  the  value  /  =  4  lies 
outside  the  instability  range. 

The  norm  for  viscous  gravity  waves  is  plotted  in  Fig.  4.  In  Fig.  4a,  A  =  0.125  as  in  Fig.  3,  and 
the  dependence  of  the  modulation  on  the  perUubation  wavenumber  is  examined.  The  amplitudes 
of  the  recurrence  are  attenuated,  as  is  also  seen  in  Fig.  3,  and  the  spikes  are  smoother.  The  value 
I  =  0.05  lies  within  the  instability  range  in  the  inviscid  limit  (A  =  0),  but  shows  decay  in  Fig.  4a. 
The  effect  of  A  is  illustrated  in  Fig.  4b  with  1  —  2  fixed.  As  A  becomes  larger,  the  amplitude 
and  period  of  recurrence  decrease;  when  A  =  0.5,  we  observe  monotonic  decay.  Fig.  4c  shows 
the  effect  of  the  disturbance  amplitude  b  when  1  =  2.  For  small  time  the  growth  rate  does  not 
demonstrate  amplitude  dependence,  whereas  for  larger  time  the  recurrence  period  increases  as 
b  decreases.  Other  calculations  for  neutral  stability  (e.g.,  1  =  4)  show  that  the  modulational 
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behavior  is  hardly  affected  by  changes  in  b  and  is  almost  exclusively  dependent  on  /,  as  expected 
from  the  linear  analysis. 

For  capillary-gravity  waves  with  surfactants,  the  dissipative  NLS  equation  is  given  by  (104). 
The  coefficient  for  the  dissipation  term  is  independent  of  the  surface  tension  T,  and  the  coefficients 
for  the  dispersion  and  nonlinear  terms  depend  only  on  T.  The  term  that  corresponds  to  the 
frequency  change  in  the  coefficient  for  the  dissipation  term  can  be  absorbed  in  the  carrier  wave 
and  does  not  contribute  to  the  magnitude  of  the  envelope  wave.  The  qualitative  behavior  of  the 
modulation  is  then  expected  to  be  identical  to  that  for  the  gravity  waves  when  T  is  outside  of 
the  range  (149),  for  the  so-called  self-focussing  type  of  NLS  equation.  When  T  is  in  the  range 
(149),  equation  (104)  is  of  the  defocussing  type.  In  this  case,  the  long-time  modulation  exhibits 
quite  different  features  from  the  simple  oscillation  with  constant  amplitude  and  frequency  found 
for  neutral  stability  of  the  self-focussing  type.  A  discussion  of  self-focussing  and  defocussing  NLS 
equations  was  given  by  Peregrine  (1983). 

Fig.  5  shows  the  evolution  of  inviscid  capillary-gravity  waves  («  =  A  =  0).  The  perturbation 
wavenumber  is  chosen  as  1  =  2,  and  four  different  values  of  surface  tension  are  considered.  For 
pure  gravity  waves  (T  =  0),  the  maximum  growth  rate  and  simple  recurrence  are  observed  as 
in  Fig.  3.  For  pure  capillary  waves  (T  =  1),  /  =  2  lies  outside  of  the  instability  region,  as  can 
be  deduced  from  (148);  thus,  an  oscillation  is  observed  with  constant  amplitude  and  frequency. 
For  T  =  0.1,  /  =  2  is  closer  to  the  lower  bound  of  the  instability  region  1  =  0,  and  so  complex 
recurrence  occurs.  The  value  T  =  0.3  lies  between  1  —  \/3/2  and  1/3,  and  the  modulation 
shows  initial  neutral  stability,  as  predicted  by  the  analysis.  However,  for  large  time  an  additional 
periodic  behavior  is  observed.  The  large-scale  period  increases  as  /  decreases,  as  can  be  seen  in 
Fig.  6a,  where  /  =  1  while  the  other  parameters  remain  unchanged  (T  =  0.3  and  k  =  A  =  0).  The 
corresponding  prcffile  for  amplitude  evolution  is  given  in  Fig.  6b.  Initially,  the  behavior  is  similar 
to  that  for  neutral  stability,  but  the  runplitude  then  gradually  decreases  and  the  frequency  starts 
to  change  until  the  initial  envelope  disintegrates  at  about  r  =  30;  still  later,  the  original  form 
of  the  emelope  is  almost  completely  recovered,  near  r  =  60.  Therefore,  this  phenomena  can 
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be  referred  to  as  recurrence,  of  a  small-amplitude  type  different  from  the  recurrence  related  to 
Benjamin-Feir  instability. 

In  Fig.  7,  we  examine  the  effect  of  dissipation  due  to  viscosity  and  a  surfactant;  A  =  0.125, 
K  =  y/2,  and  /  =  l/y/6.  The  perturbation  wavenumber  corresponds  to  the  maximum  growth  rate 
y/ff  for  in  viscid  pure  capillary  waves.  The  norm  for  the  case  T  =  1  grows  to  a  maximum  at  twice 
its  initial  value  near  r  =  50  (Fig.  7a).  When  T  =  0.1,  the  instability  with  complex  recurrence  is 
observed  despite  viscous  dissipation,  whereas  for  T  slightly  smaller  than  1  —  y/3/2  (T  =  0.1339), 
the  complex  recurrence  is  suppressed  by  dissipation,  and  monotonic  decay  is  observed.  For 
T  =  0.3,  a  disintegration  of  the  wave  envelope  is  obvious  even  with  viscous  damping,  as  can  be 
seen  clearly  in  the  corresponding  evolution  profile  in  Fig.  7b.  FVom  Fig.  7,  we  can  also  deduce 
that  the  period  of  the  disintegration-recovery  has  increased  from  that  in  Fig.  6  because  a  smaller 
i  has  been  considered. 

In  all  the  calculations  above,  the  initial  condition  used  is  (152),  which  corresponds  to  a 
symmetric  amplitude  modulation.  We  now  consider  more  general  modulations  by  using  (121) 
instead,  with  Bi(0)  and  ^3(0)  not  necessarily  equal  and  real  as  in  (152).  For  simplicity,  the 
results  for  inviscid  gravity  waves  are  presented  in  Fig.  8.  Only  the  dominant  fundamental  and 
the  sideband  modes  are  plotted.  In  all  the  cases  considered,  the  wavenumber  perturbation  is 
/  =  2,  so  that  all  exhibit  an  initial  Beqiamin-Feir  instability  followed  by  simple  recurrence.  In 
Fig.  8a,  the  given  initial  condition  is  identical  to  that  in  Fig.  la  and  b,  and  so  the  results  are 
identical.  In  Fig.  8b,  it  is  obvious  that  the  phase  difference  changes  the  modulational  behavior 
significantly,  including  the  recurrence  period.  Fig.  8c  shows  the  effect  of  different  initial  sideband 
amplitudes,  here  showing  an  increase  in  the  recurrence  period  with  a  decrease  in  amplitude.  In 
Fig.  8c,  only  one  of  the  sideband  inodes  is  present  initially,  but  subsequent  evolution  shows  that 
the  other  mode  is  automaticaUy  excited  to  produce  the  BenJamin-Feir  instability.  The  initially 
different  sideband  amplitudes  in  Fig.  8c  become  nearly  identical,  as  predicted  by  Stiassnie  ii 
Kroszynski  (1982),  but  they  recover  their  difference  periodically  in  the  subsequent  recurrence. 
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7  Concluding  remarks 


The  methods  of  multiple  scales  and  matched  asymptotic  expansions  have  been  used  in  a 
formal  derivation  of  evolution  equations  for  weakly  nonlinear  viscous  water  waves.  The  result  is 
the  nonlinear  Schrodinger  equation  with  an  addition^  linear  dissipation  term.  For  both  gravity 
waves  and  capillary-gravity  waves  with  an  insoluble  surfactant,  the  largest  terms  in  the  damping 
coefficients  are  identical  to  classical  linear  results.  For  capillary-gravity  waves,  nonuniformities 
are  observed  despite  dissipation  due  to  the  surfactant.  Near  the  second-harmonic  resonance, 
a  modified  set  of  evolution  equations  is  obtained,  which  are  matched  asymptotically  with  the 
nonresonant  results.  For  pure  capillary  waves,  a  rescaling  of  the  long-scale  free-surface  elevation 
is  required  in  the  presence  of  curfactants. 

The  derived  evolution  equations  are  solved  analytically  and  computationally  to  examine  the 
modulation  of  infinite  wavetrains  with  sideband  disturbances  of  the  Benjamin-Feir  type.  The 
linear  analysis  shows  that  in  the  presence  of  dissipation  the  modulation,  described  by  modified 
Bessel  functions  of  complex  order,  has  the  same  initial  behavior  as  the  inviscid  case  but  eventually 
decays  to  zero.  In  the  inviscid  limit,  the  instability  condition  obtained  originally  by  Beiyamin 
&  Feir  (1967)  is  recovered  for  gravity  waves,  whereas  a  stable  range  of  Weber  number  is  found 
for  capillary-gravity  waves,  which  agrees  with  that  obtained  by  Djordjevic  k  Redekopp  (1977). 
In  this  range,  spectral  computations  show  a  new  type  of  recurrence  not  directly  related  to  the 
Benjamin-Feir  instability. 
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LIST  OF  FIGURES 

Figure  1.  Evolution  of  modulations  for  inviscid  gravity  waves  (A  =  0)  when  the  ampUtude  per¬ 
turbation  6  is  0.1:  a)  /  =  1,  b)  /  =  1  Fourier  space,  c)  /  =  2,  d)  /  =  2  Fourier  space,  e)  /  =  4,  f) 
/  =  4  Fourier  space. 

Figure  2.  Evolution  of  modulations  for  gravity  waves  expressed  by  N{t)/N(0)  for  /  =  1 
( - ),  /  =  2  ( - ),  /  =  2n/2  ( . ),  and  /  =  4  ( - ). 

Figure  3.  Evolution  of  modulations  for  gravity  waves  with  small  viscosity  (A  =  0.125):  a) 
/=l,b)/  =  2,  c)/  =  4. 

Figure  4.  Evolution  of  modulations  for  gravity  waves  when;  a)  A  =  0.125  amd  b  =  0.1  are 

fixed,  and  the  perturbation  wavenumber  varies  for  /  =  0.05  ( . ),/=!( - ),  f  =  2  ( - ), 

and  /  =  4  (—  •  — )  b)  /  =  2  amd  6  =  0.1  are  fixed,  and  the  viscous  dissipation  varies  for  A  =  0 

( - ),  A  =  0.125  ( - ),  A  =  0.25  ( - ),  A  =  0.5  ( . ),  and  A  =  1  ( - )  c)  /  =  2  and 

A  =  0.125  are  fixed,  and  the  amplitude  perturbation  varies  for  b  =  0.01  ( - ),  6  =  0.1  (-  -  -  -), 

6  =  0.2  ( - ),  and  6  =  0.3  ( . ). 

Figure  5.  Evolution  of  modulations  for  inviscid  ci4>illary-gravity  waves  (A  =  k 
the  perturbation  wavenumber  1  =  2  and  the  Weber  number  is  chosen  as  T  =  0  ( — 

(----),T  =  0.3  ( - ),  andr=:l( . ). 

Figure  6.  Evolution  of  modulations  for  inviscid  capiUary-gravity  waves  (A  =  k  =  0)  when 
the  perturbation  wavenumber  /  is  1  and  the  Weber  number  T  is  0.3  expressed  by:  a)  N(r)/JV(0) 
b)IA((,r)l. 

Figure  7.  Evolution  of  modulations  for  c^illsry-gravity  waves  in  the  presence  of  viscosity 


=  0)  when 

-).r  =  o.i 
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(A  =  0.125)  and  8urf2u:tant  (k  =  y/2)  when  the  perturbation  wavenumber  /  =  l/y/S  and  the 

Weber  number  is  chosen  as:  a)  T  =  0.1  ( - ),  T  =  0.1339  T  =  0.3  ( —  ■  — ),  and 

T  =  1  ( . )  b)  r  =  0.3. 


Figure  8.  Evolution  of  dominant  Fourier  modes  (n  =  0  and  ±1)  for  inviscid  gravity  waves 
when  1  =  2.  The  initial  sideband  amplitudes  are  represented  by;  a)  Bi  =  0.05  and  =  0.05,  b) 
Bi  =  0.05  and  B2  =  0.05 1,  c)  Bi  =  0.05  and  B2  =  0. 
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